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DESIGN PROCEDURE AND LIMITED TEST RESULTS FOR A HIGH SOLIDITY, 12 -INCH 
TRANSONIC IMPELLER WITH AXIAL DISCHARGE 
By Linwood C. Wright and Karl Kovach 


SUMMARY 

An approximate inverse impeller design procedure based on the 
assumption of axisymmetric, nonviscous flow is presented along with some 
l imi ted experimental data obtained from a high solidity, 12 -inch tran- 
sonic Impeller. The design procedure consists of a point -by -point solu- 
tion of the Impeller flow conditions along a series of streamlines in 
the radial-axial plane. Also, a method is presented for estimating the 
blade loading from the moment of momentum considerations when the mean 
channel flow is assumed parallel to the mean blade surface and the static 
pressure is assumed to vary linearly across the passage. 

The transonic impeller designed for an equivalent tip speed of 
1200 feet per second by means of the subject procedure consists of a 
0*40 hub-tip ratio, swept -back leading edge, aluminum-b laded rotor. 

The design performance was not obtained experimentally because of a 
restriction of the flow in the outlet annulus downstream of the impeller 
which causes the flow to be about 14 percent less than the design value. 
The resulting large angles of attack at the blade leading edges caused 
vibrations which induced partial structural failure before alteration of 
the outlet passage could be accomplished. After the initial partial 
fracture of one of the blades, total failure was delayed by installation 
of damping wires threaded through the blades. Although further running 
time could be obtained by use of damping wires, the life of this config- 
uration was too short to justify extensive alteration of the outlet in 
order to alleviate the choking. 

Because the flow was choked in the outlet annulus, at design equiv- 
alent tip speed the experimental performance was restricted to a pressure 
ratio of 2.164 at an equivalent weight flow of 26.47 pounds per second 
(33.71 lb/sec-sq ft of frontal area) and an impeller adiabatic effi- 
ciency of 81 percent . 

The impeller was designed to have a static-pressure drop from inlet 
to exit in order to reduce the adverse viscous effects. However, the 
losses due to resulting high stator inlet Mach numbers are concluded to 
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outweigh any favorable effects on the boundary layer. Also, the advis- 
ability of approaching equal energy input along the blade span is 
indicated. 


INTRODUCTION 

For some time compressor designers were of the opinion that the 
performance of axial-flow compressors would be limited by a maximum 
blade relative Mach number. The belief existed that after the fashion 
of isolated subsonic airfoils, performance deteriorated rapidly once 
the blade relative speeds approached the critical or force-break values 
(ref. l) . Experiments indicated that this belief was reasonably accurate 
for the conventional airfoil sections and rotor configurations utilized 
in subsonic axial-flow compressors. 

One -dimensional analysis of supersonic flow configurations indicated 
that if the extended bow wave pattern observed when subsonic blading was 
placed in a supersonic stream could be eliminated and the normal shock 
confined within the passage, then reasonable efficiencies and consider- 
ably increased pressure ratios might be realized. The proposed method 
for eliminating the extended wave pattern incorporated the earlier prin- 
ciple of the Busemann biplane to contain the wave system wholly within 
the blade passage. This requirement led to the necessity of maintaining 
the entrance relative Mach number sufficiently large to allow passage of 
all of the flow for the given blade contraction without spillage (refs. 2 
and 3) . 

Early supersonic compressors were designed with the idea of meticu- 
lously avoiding blade relative hub Mach numbers below those which would 
permit shock entrance at design speed along the entire blade span. This 
restriction placed a limit on the minimum hub -tip radius ratio for a 
given inlet velocity and rotational speed and hence limited the maximum 
weight flow. This minimum Mach number restriction led to the use of 
relatively high tip speeds (usually about 1500 or 1600 ft/sec) and the 
attendant severe structural problem. Moreover, attention was focused 
away from the intermediate tip speed range (1100 to 1400 ft/sec) where 
reasonable radius ratios and resulting large weight flows were struc- 
turally feasible. This speed range was avoided because the resultant 
relative Mach numbers fell in the transonic speed range where the extended 
wave pattern could be neither eliminated nor analyzed. 

Experimental results of several supersonic compressors, however, 
indicated that their performance in the speed range described as tran- 
sonic (1100 to 1400 ft/sec) was in general better than the design speed 
performance (see refs. 4 and 5) • The discontinuous variation in flow 
conditions with speed which might have been expected (at the speed of 
shock entry) was not observed. There was evidence that large losses 
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were not inherent in working with transonic flows. In the absence of 
viscous wall effects this might have been readily predicted inasmuch as 
the shock losses are inherently small. 

The above observations led to consideration of the transonic com- 
pressor configuration which is characterized by subsonic blade inlet 
relative velocities at the hub and supersonic blade inlet relative veloc- 
ities at the tip at rated speed. The configuration is such that arbi- 
trarily low hub -tip radius ratios may be utilized with greatly increased 
weight flow per unit frontal area. The theoretical pressure ratio per 
stage generally can exceed considerably that of a subsonic stage while 
falling somewhat below that of a supersonic stage. The efficiency will 
theoretically exceed that of a supersonic stage but probably cannot equal 
that of a subsonic stage of optimum design. The generally reduced tip 
speed allows the use of increased blade span for a given stress level. 

Thus, an analytical and experimental investigation of the compres- 
sor type termed transonic was undertaken for the purpose of developing 
an appropriate design procedure, designing, testing, and analyzing the 
performance of a high solidity transonic impeller. This impeller differs 
from the axial-flow transonic impeller of reference 6 primarily in that 
the current impeller utilizes a lower hub-tip ratio, more turning, a 
considerably higher solidity, and a decreasing mean static pressure 
through the impeller. The accelerating mean flow was intended to pre- 
vent large boundary layer growths and eliminate the necessity for large 
boundary layer allowances to prevent choking. A necessary consequence, 
however, is the undesirably high Mach number level of the impeller 
absolute leaving flow. 

This report presents the transonic impeller design procedure and 
limited test results. The experimental results are incomplete because 
of inability of the test rig to pass the design weight flow and are 
insufficient to warrant strong conclusions regarding the questionable 
advantages of an accelerating mean flow through the rotor. The occur- 
rence of structural failure did not allow sufficient time to make the 
changes in the exit flow passage which would have been required to 
permit passage of the design flow. 


DESIGN PROCEDURE 

The first objective of the program was the development of a design 
procedure. Complete freedom in the aerodynamic specifications would be 
highly desirable; however, for the desired configuration of relatively 
high wheel speed and long blade height, structural considerations will 
necessarily restrict the aerodynamic freedom. 
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In order to control centrifugal stresses, specification of the 
radial distribution of blade thickness should be allowed. In addition, 
radial blade elements will be desirable to eliminate excessive bending 
stresses. The only remaining freedom that is available to satisfy aero- 
dynamic specifications is then the blade mean line angle and thickness 
distribution along a single streamline or specified contour in the 
radial-axial plane. The design procedure under these conditions reduces 
to two parts: (l) The selection of the blade shape and thickness dis- 

tribution along a given streamline, and (2) the determination of the 
rotor hub and tip annular boundaries that permits the required weight 
flow to pass through the impeller compatible with the desired conditions 
existing on the specified streamline. The remaining portion of the 
blade other than the specified section will be fixed by the condition 
of radial blade elements and a specified radial taper. 

The procedure herein described is usually designated as the inverse 
procedure in compressor design practice. That is, the tip speed and 
performance are specified and the necessary impeller geometry is deter- 
mined. 


Simplifying Assumptions and Discussion 

The following are the simplifying assumptions: 

1. The flow processes throughout the impeller are considered isen- 
tropic (a displacement thickness allowance for boundary layer growth may 
be made, however) . 

2. The flow is assumed axially symmetric and the mean flow path is 
assumed parallel to the blade mean lines. 

3. The continuous flow variations in the impeller are presumed to 
be adequately represented by incremental changes along the radius r 
and the axis z. (See list of symbols in the appendix.) 

That the entire flow will not be isentropic is apparent. If, how- 
ever, the pressure gradients can be maintained below the limit which will 
cause separation and boundary layer displacement thickness allowances are 
made thereby preventing choking and distorting of local stream tubes, the 
free-stream flow paths should not differ essentially from those isen- 
tropically computed. In the absence of large adverse pressure gradients 
the methods of reference 7 might be conveniently utilized to calculate 
the approximate boundary layer growth along any surface where the free- 
stream conditions are known. 

The commonly used assumption of axial symmetry is strictly true 
only for an infinite number of zero thickness blades. This condition is 


MCA RM E53B09 


5 


believed to be reasonably approximated for high solidity thin blades such 
as those of the present investigation (see ref. 8). Care should be 
exercised, however, and possibly special allowances made in the continu- 
ity equation when the blade mean relative Mach number approaches 1. 

Under these conditions the use of a sonic mean velocity in the continu- 
ity equation will overestimate the flow capacity of a passage with mean 
velocity in the transonic range. 


Derivation of Equations 

Nomenclature and coordinate system . - Throughout the current pro- 
cedure cylindrical coordinates were used with the nomenclature illus- 
trated in the following sketch. The relations between the various veloc- 
ity components were developed as follows : 


+z 



COS P ' = V^V' g 

(1) 

sin cp / = V r /V’ 

(2) 

tan cp = V r /V z 

(3) 

( v ) 2 = ( v ^ 0 ) 2 + ( v r ) 2 

(4) 


From equations (l) and (4), 



( 5 ) 
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From expressions (5) and (3) , 

v 2 _ (V cos 3') 2 ( 6 ) 

v z “ 2 2 

1 + tan cp cos P 1 

Determination of blade mean camber line . - Utilization of an ana- 
lytical expression for the blade mean surface in three dimensions is 
convenient in the computations of the impeller flow. The following 
development with the aid of the sketch below relates the radius, axial 
position, and the angles p ' and 9 as required for radial blade 
elements: 

. . A 6 d 6 tan P ' 

lim - 7 — = "3 — = — 

. Az dz r 

Az — *0 



All linear dimensions may now be made dimensionless by dividing by 
the tip radius r-f Then d0/d(z/rt) may be equated to a polynomial in 
z/r-fc, such as in the following expression: 


d0 

d(z/r t ) 



tan p 1 
r/rt 


(7) 


In order to determine the coefficients a, b, and c, three condi- 
tions must be specified. The following three conditions are proposed 
for convenience : 
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1. At z/r-t = 0, C = tan 3 ’/r/r-fc where 3 ' may be found as 

3' = tan - - 1 - oar/v^-t- for no guide vanes. The value of V z , which is the 

axial velocity component at the tip, is known from the inlet computations. 

2. At (z/r^-Jg, 3' may be specified at a particular radius, for 
example, 3^ = 0. 

3. The third condition may be selected so as to distribute the 
blade loading relatively uniformly over the blade surface. In the pre- 
sent design, a fixed angle corresponding to a predetermined degree of 
turning was specified for z/r^ = 0-30. In this manner, the peak load- 
ing was placed near the leading edge for the blade defined by the coef- 
ficients specified. 

Blade thickness distribution . - A blade thickness distribution is 
specified along an initial streamline. For structural reasons a blade 
radial thickness taper with the maximum thickness at the hub is desir- 
able. With the mean surface equation known, if a particular radial 
taper is set, the specification of the thickness distribution along one 
streamline fixes the entire blade geometry except for the unspecified 
radial boundaries. 

Determination of pressure distribution around a blade element . - 
In order to evaluate the aerodynamic desirability of any given blade 
shape as determined by a mean line and a thickness distribution obtained 
from the preceding procedure, an estimate of the static -pressure or 
blade loading distribution is required. An approximate method for 
determining the static -pressure distribution follows. This procedure 
is a variation of that used in reference 9 utilizing moment of momen- 
tum principles. Though basically identical, the methods differ in two 
particulars. The current procedure involves the determination of 
blade surface velocities for a predetermined blade shape and the 
assumption of a linear variation of pressure between blade driving- 
and trailing-face surfaces. In reference 9 the blade loading diagrams 
were selected as desired and the blade shape was evolved from moment 
of momentum considerations; the velocity variations between blades were 
assumed linear. 

The static -pressure difference across a blade passage in the tan- 
gential direction may be approximated by equating the rotor torque to 
the rate of change of angular momentum of the fluid. Then for con- 
stant radius, 


torque = Mass • r -A Vg 

where AVg is the change in the absolute tangential velocity over the 
surface interval along the streamline for which the torque is computed. 
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Now 


As 


n 


Ap q X r X h X dz 


vj 


s n-l 


2jtr-j_ 

— H — h iPi v i cos Pi cos ^ AV 0 X r 


( 8 ) 


where AVg is the change in Vg between stations n-1 and n in the 
radial-axial plane and this increment is selected sufficiently small 
that r may he considered constant, and Ap 0 is the net pressure act- 
ing in the tangential direction on the blade surface of height h 
between the stations n-1 and n. 


In order to compute Ap q at any point z, the expression (8) may 
be written in the finite difference form below assuming a linear varia- 

tion in the quantities for small differences The linear dimen- 

sions are made dimensionless through division by the tip radius r^ and 
aor is in feet per second. 



(9) 


-) 

r t/n- 


(V z0 sin p« - cnr) 


n-1 


This expression may be used within the accuracy of the current 
assumptions by utilizing the radial average of the upper and lower stream 
surface conditions along the mean radial position of the incremental 
annular segment. As a consequence of the assumption of axial symmetry, 
stream conditions are further averaged between the blade driving- and 
trailing-face values. In the above expression a single bar refers to the 
average or mean value between the appropriate boundaries, radial or tan- 
gential, and a double bar refers to conditions averaged in both 
directions . 

The problem of obtaining the exact loading from the pressure dif- 
ference Ap q has not yet been adequately solved. There are several 
assumptions, however, by which solutions may be obtained. For the cur- 
rent design a linear variation of pressure from blade driving to trailing 
face is assumed. The following expressions result for surface pressures: 

Pdf = P + Ap e /2 
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p tf = p - Ap e/ 2 


The blade driving- and trailing-face relative velocities may be 
readily obtained from the preceding surface static pressures in a manner 
which will be described later. 


Determination of radial pressure gradient . - Within the assumption 
of axially symmetric, nonviscous flow the radial forces acting on the 
fluid particles passing through an impeller of the type where radial 
velocity components are relatively small may be approximated as shown 
below. 


In references 10 and 11 (pp. 990-997) the radial component of 
Lorenz ' s equation of motion is given as 

F _l*P + d!r + V > 2 

r p 3r d t 2 r 

! dp (“*• + V) f 

However, for radial blades F r = 0 and — ^ = 


d 2 r 

dt 2 ‘ 



From the preceding sketch, which shows the projection of a stream- 
line in the radial -axial plane, dcp = — ^ or ^ Pn Edition, 

v = or dt = where 's' is the distance along a streamline 

projection in the radial-axial plane. 
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Now 


or 


1 dp 

p dr ” 


1 dp 

p ¥" 


^cor 


9 


(v rz cos cp ^ + sin cp 


(cor + V e ') 2 V ; 


rz 

R„ 


dV 


dt / 


rz 

cos cp - v rz Sin cp 


Making all lengths dimensionless gives 

>2 „ 2 


dp 


Pa £) 


(cur + V 0 ')‘ 
r 

r t 


rz 

Rc 

r t 


dV. 


cos cp - V. 


rz 




sin qp 


( 10 ) 


The expression for radius of curvature of a smooth curve is used to 
evaluate 
plane . 


R c , the streamline radius of curvature in the radial-axial 


R „ - 


[* + (ml 


3/2 


d 2 r 

dz 2 


( 11 ) 


In order to utilize equations (10) and (ll) for obtaining the radial 
pressure gradients } reference is made to the third assumption that con- 
tinuous variations are assumed to be adequately approximated by small 
finite incremental changes. Thus dr/dz = tan cp and 

2 /,_2 tan - tan %-l 


d 2 r/dz^ s 


Az 


where subscripts 


n-1 and n are succes- 
will have the sign 


sive points along the streamline \|r - constant. 

d^r/dz^ and will become the dimensionless term R c / r t w k en Az is 

2/2 

made dimensionless through division by r-t in the expression d r/dz . 
Again considering the change along \|r = constant from n-1 to n, 


dV. 


rz 


AV 


rz 


(^rz^n “ ^rz^n-1 


‘ <*> £). ■ (*),., 


where s/r-^ is the dimensionless distance along the stream surface 
projection in the radial-axial plane. 
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Computation of radial boundary from continuity . - The continuity 
relation to be satisfied everywhere between the stream surfaces is 


V 1 “ P]_V]_ cos ^1 ^ 


2rtrj 

"N - 


hi = 


p y ' cos v cos p' (rr - h 


from which 


_h_ 

^1 


p l V l 


cos cp 



pV COS cp' COS P ' 


( 


2n JL 
r t 
N 



( 12 ) 


where subscript (l) refers to the rotor inlet and the double bar indi- 
cates the mean tangential and radial positions. The linear dimensions 
are made dimensionless on division by radius r-^. Other terms .are 
defined in the list of symbols. The angles are approximated as follows: 




A convenient form of the energy equation as developed in reference (12) 
is given and may be used for computations of the stream density: 


pt,o I 6 




(13) 


where Gj q is the stagnation velocity of sound in front of the 
impeller . 


Determination of velocities from energy relations . - In order to 
obtain the mean velocity distribution on a new streamline from the known 
velocities on an adjacent streamline and the known Ap r distribution 
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between the streamlines, the energy equation and the adiabatic relation 
may be utilized. The energy equation is given in a convenient form for 
no prewhirl. 


( v ' ) 2 

K m+1 7 . n + 


- _ 2 H 

(V 1 ) 

m i r t 

Kn+i )2 - k/ 

2 + 


2 + ^p in 

2 


or 

(=Vi ) 2 = (a* m+1 ) - (^ m ) 2 + 2C ( t m - + (V 'm ) 2 

where the subscripts m and m+1 refer to the adjacent streamlines 
\|r = m \|r = m+1. The double bar indicates the mean radial position 

and the mean tangential position between the driving and trailing faces. 

If " ^m+l “ (At) r > 

(V m+ p) = (^m+p) ~ (^m ) + 2Cp(At) r + (V' m ) (14=) 

For the case where the computations proceed inward toward the axis, 

(Atl = t , . - tn If the flow is considered 

u upper mean line lower mean line 

isentropic the temperature at m+1 may be found from (Ap) r using the 
adiabatic expression 


r-i r-i 



The right side of this equation may be replaced by the first three terms 
of a binomial series expansion with essentially no error. 

General computational procedure - - After a satisfactory blade mean 
line and thickness distribution have been established, the hub and tip 
contours required to pass the desired weight flow and satisfy the con- 
ditions prescribed along the initially specified streamline must be 
obtained. The proposed procedure may be considered as a stream filament 
method applied through the inlet section and the impeller . The flow at 
the entrance to the inlet section is first divided into several segments 


NACA RM E53B09 


13 


of radial height . With the specified streamline as one boundary, 
the variation in h for the adjacent stream tube through the entrance 
section and rotor (along z) necessary to satisfy continuity is obtained 
from equation (12). The radial gradient in pressure and the velocity on 
the subsequent streamline are next computed from the flow on the speci- 
fied streamline with the aid of equations (10) and (14) . The pressures 
along the adjacent streamline are estimated by extrapolating the pressure 
gradients over the distances h obtained from continuity. A more accu- 
rate estimate of the mean annular conditions may now be obtained from 
the values on the upper and lower boundaries. These mean conditions may 
now be used to obtain the final orientation of and the conditions along 
the second streamline. This process may be continued progressing from 
the specified streamline in either or both radial directions until the 
desired mass flow is contained. 

The blade loading may be estimated from expression (9) derived from 
moment of momentum considerations. 


The expressions previously derived or presented may be utilized in 
the detailed computation of the inlet and impeller configuration in the 
following manner. 

Entrance section . - The entrance conditions may be computed by means 
of the same relations utilized in the rotor. The rotor equations reduced 
to the absolute reference frame for use in the inlet are as follows: 

From continuity 


Upon elimination of blade considerations the radial pressure gra- 
dient becomes 


Application of Equations to Compressor Design 



where 


1 
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1 dp 



V rz 2 cos cp 
Rc 


r t 



sin cp 


The energy equation yields Mach number and hence velocity directly 
in the following form: 



In a manner identical to that described under the general procedure, 
the velocities may be specified on a single streamline in the entrance 
region starting from some position upstream of the impeller and ending 
at the impeller blade leading edge. By computing the streamlines radi- 
ally away from the one specified, as was described for the impeller, 
the i nnpr and outer walls upstream of the impeller may be computed. 


Matching of entrance and blade leading-edge conditions . - If the 
rotor blades are restricted to radial elements, as is done in equa- 
tion (7), the blade angle will match the inlet flow angle at a constant 
axial position (z » 0) at all radii only for the case of a constant 
entrance velocity with no prewhirl. 

If an inflow velocity other than constant is desired at the leading 
edge of a blade with radial elements, a swept leading edge is required 
to match the inflow and blade angles at all radii. The following semi- 
graphical procedure may be utilized in obtaining the blade leading-edge 
contour : 


In the region near the leading edge, the relative flow angle, 
obtained from the rotational velocity and the velocity determined from 
the solution of the inlet flow along the streamline, is computed for 
several stations along z. From the values of r/r-^ obtained for the 
streamline in the inlet section, the corresponding blade angle is com- 
puted, using the prescribed equation for the blade mean surface. The 
z-position at which the relative inflow angle matches the blade angle 
establishes the blade leading edge for zero incidence angle. 


Blade leading-edge thickness. - The leading-edge thickness may be 
specified arbitrarily or from structural requirements. Usually, in 
order to minimize the leading-edge disturbance, the smallest practical 
leading-edge angle and thickness are utilized. The leading-edge angle 


2 

may, in general, be defined as 2 tan ^ where T is the blade 

thickness on the subject streamline and is located A(z/r-t) behind the 
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blade leading edge. At this point T may be determined from the radial 
thickness taper and the blade thickness at the corresponding z-position 
on the previous streamline. 

When this procedure is applied to compressors having subsonic or 
transonic relative inlet velocities some deviation from smooth inlet 
flow is probably unavoidable, with the magnitude of the disturbance a 
function of the perturbations caused by the blades. For supersonic 
blades, however, the disturbances may, if desired, be contained within 
the blades by suitable selection of the number of blades and their 
thickness . 


COMPUTATIONAL EXAMPLE AND DESIGN RESULTS 

In order to illustrate the order of the computational procedure, 
the following abbreviated example of the process followed in the design 
of the subject impeller is presented along with the design results. 

The following constants were initially specified: 

D^ = 12.00 in. 

U t = 1200 ft/sec 

P' 2 = 0 for all radii 

Djj/D t = 0.40 at the inlet 

Axial depth = 6-00 in. at the hub 

M-fc = 0.97 at the inlet 

Blade mean line specification . - A cylindrical outer shroud was 
selected and the initial streamline was specified on this surface. The 
blade mean surface, which under the assumption of axial symmetry is 
identical with the mean relative stream surface, is obtained from equa- 
tion (8) where the constant coefficients are fixed by the following 
conditions on the tip streamline. 

1- p 2 . = tan - ' 1 ' cur t /V z 

2 . P 2 =0 

3. p* (z/r t = 0.30) = 35° 

The blade mean surface is thus represented by the equation 

o 

2.354 (z/r^.) - 4.060 z/r^. + 1.706 - tan P'/r/r^ 
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Specification of "blade thickness- - The blade thickness distribu- 
tion about the tip mean line previously defined was arbitrarily fixed at 
the minimum values believed to be structurally sound. The blade thick- 
ness measured normal to the blade mean line is given in figure 1- 

The thickness taper was defined by a blade included angle of 1° 
between the driving and trailing face blade elements in the plane normal 
to the blade mean line at that point. The blade thickness at any radius 
is then obtained from the corresponding blade tip thickness as 

T/rt = (t/rt)t + [( r / r t)t - (r/rt)] 0.01714 

No allowance for boundary layer was made in these calculations. 
However, a displacement thickness allowance may be conveniently applied 
to either the thickness distribution or the computed value of h if 
desired. 

Computation of radial boundaries . - Equations (ll) and (12) were 
utilized to compute the annular segmental boundaries progressively from 
the tip cylindrical boundary to the hub contour- The resulting hub 
contour gave a hub radius ratio change from 0.40 at the inlet to 0.619 
at the exit. Equations (13), (14), (15), and (16) were used to obtain 
the radial pressure gradients and velocity as these quantities are 
necessary to the boundary computations. 

Determination of blade sweepback . - In the process of computing 
the radial boundaries, the location of the impeller blade leading edge 
is accomplished in the manner previously described by matching the blade 
mean surface to the inlet relative flow direction. The total blade 
leading-edge sweepback of the subject impeller is found to be 

(z/r t ) t - ( z / r t)H = °' 197 

Computation of blade loading characteristics . - The radial average 
of the blade loading in an annular segment of height h may be .esti- 
mated from a solution of equation (10). The tip loading in terms of the 
driving and trailing face velocities is plotted in figure 2. 

Design performance. - The design resulted in the following theo- 
retical p'erf ormance : 

Over-all pressure ratio of impeller less normal shock losses = 2.41. 

Maximum design speed weight flow = 30-8 lb/sec (39.2 lb/sec-sq ft of 
frontal area). 

Over-all adiabatic efficiency considering only normal shock losses = 
0.852. 
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The tip inlet velocity diagram is shown in figure 3, while fig- 
ure 2 gives the comparable loading diagram and mean velocity distribu- 
tion. The restriction to a small rise in mean velocity through the 
impeller was intended to reduce the amount of deceleration on the blade 
trailing face at the critical trailing-edge portion. In spite of this 
unorthodox mean pressure decrease through the impeller, the trailing - 
face velocity (fig. 2) still undergoes a substantial deceleration as the 
flow approaches the trailing edge. The effects of this deceleration on 
wake thickness and separation tendencies are not known. 

Actually only the loading trends can be obtained from the moment of 
momentum; thus, it does not appear likely that the blade would unload 
before the trailing edge as indicated in figure 2. This phenomenon 
results from an inflection point in the blade mean line to which the 
relative rotor flow is assumed parallel. 


EXPERIMENTAL APPARATUS AND TEST PROCEDURE 

The transonic impeller which was machined from a solid 14S-T alu- 
minum forging is pictured in figure 4. It was tested in the 3000 horse- 
power rig shown schematically in figure 5. The piping provided for 
refrigerated or atmospheric inlet air as well as altitude exhaust. A 
calibrated submerged adjustable orifice was used to measure the weight 
flow. Most of the tests were run at 15 inches of mercury absolute inlet 
pressure. No differentiation is made between these tests and those few 
made with atmospheric inlet pressure, since no distinct differences could 
be observed in the reduced data. 

The instrumentation consisted of temperature and pressure probes 
(fig. 5) plus static-pressure orifices located along the outer casing 
before, over, and after the impeller. Static-pressure orifices were 
also located along the inner wall (hub) before and after the impeller. 
Five calibrated double stagnation thermocouple probes were located 
6 inches downstream of the impeller trailing edge. The annulus area was 
divided into five equal annular rings and a probe located in the center 
of each ring. Six spike-type thermocouples, six total-pressure, and two 
static -pres sure probes were located in the inlet depression tank where 
the total area was such that the static and total pressures were 
essentially eqLual. 

A probe actuator was located downstream of the impeller as indi- 
cated in figure 5 so that either a combination cone probe (static and 
total pressure plus yaw) or a claw total -pressure probe could be used 
approximately 1 inch behind the impeller blade trailing edge. 
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Both probes were calibrated, for angle in a subsonic and. supersonic 
stream. The cone probe was calibrated for angle and static pressure. 
For Mach numbers above approximately 1.05, a close agreement was noted 
between the computed surface static pressure and that measured on the 
surface of the cone. When the claw total probe was used, the stream 
static pressure was obtained by interpolation between the inner and 
outer wall orifice measurements- Ten survey points were usually taken 
across the approximately 2 inch impeller exit annulus. 

Accuracy of data. - Observed measurements are estimated to be 
accurate within the following limits : 


Temperature, °F . . . • 
Pressure, in. Hg • • • 
Air weight flow, percent 
Impeller speed, percent 


± 0.8 

±0.05 

± 1.0 

±0.5 


On the basis of these tolerances, the maximum errors in the per 
formance figures are approximately 


Pressure ratio (at lowest value) 

Adiabatic efficiency (percent at lowest temperature rise) • • • 


No tolerances are given for the angular measurements inasmuch as 
they have small effect on the total-pressure ratio and adiabatic effi- 
ciency and moreover are difficult to estimate with any degree of assur- 
ance in the impeller wakes. 

In the table of performance data (table i) two adiabatic efficiency 
points at a tip speed of 500 feet per second have values equal to or 
slightly exceeding 1-00 (1.00 and 1-03). The combination of temperature 
and pressure errors giving the maximum adiabatic efficiency error would 
reduce the recorded efficiency only 2-5 percent. There is one other 
source of error, however, which is not exactly a function of the accuracy 
of measurements. During the testing period the necessity of holding con- 
ditions constant for a considerable period was found necessary in order 
to obtain stable temperature readings. This phenomenon has been commonly 
noted where there is high heat capacity. Hence the theory is advanced 
that those efficiency points exceeding 1.00 have an additional small 
error because of insufficient time allowance for the measured temperature 
to reach its maximum value for the fixed conditions. 


Within the limitations specified, the data provided by this instru- 
mentation were reduced to the performance data presented. 
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EXPERIMENTAL RESULTS AND DISCUSSION 
Over -all Performance 

Presentation of the 12 -inch transonic mixed-flow impeller perform- 
ance characteristics is complicated hy a data deficit resulting from 
inability to obtain the high weight flow side of the performance map. 
Because of these conditions care was exercised in an attempt to elimi- 
nate erroneous conclusions which might result from the inadequate data. 

The measured over-all performance is given in table I and in figure 6 which 
shows adiabatic efficiency contours superimposed on a plot of impeller 
pressure ratio against equivalent weight flow for 500, 700, 800, 1000, 

1100, and 1200 feet per second equivalent tip speeds. It is apparent 
that only the low weight flow portion of the curve is shown in figure 6. 

The shape of the curves indicates that the peak impeller pressure ratio 
was not attained. That part of the performance map determined at all 
speeds above 800 feet per second was considerably below design expecta- 
tions. There are several indications, which will be discussed later, 
that choking in the annulus after the impeller is responsible for the 
failure to obtain the weight flow corresponding to peak pressure ratio. 

Observation of the pressure ratio weight flow curves (fig. 6) shows 
that only the lowest tip speed (500 ft/sec) indicates pressure ratio and 
efficiency peaks at less than the maximum measured flow; all other curves 
display continually rising pressure ratio and efficiency with weight flow 
up to a maximum. Moreover, there is a relatively small increase in weight 
flow with speed which is reflected in the large increase in minimum angle 
of incidence with speed shown in figure 7. It might be expected that, 
for minim um back pressure on the rotor, the axial inlet velocity (hence, 
weight flow) would increase with speed to an extent at least sufficient 
to maintain nearly constant impeller blade incidence angle. This char- 
acteristic, which should be exhibited even by a choked impeller, is not 
apparent. Instead, the incidence angle (fig. 7) increases with speed. 

This condition is clearly not characteristic of an unrestricted flow 
path. 


Impeller Flow Limits 

Impeller inlet flow . - Utilizing the current design procedure, the 
resulting transonic impeller blade leading edges were swept behind the 
hypothetical Mach cone based on the theoretical relative inlet Mach num- 
ber. The structural and aerodynamic effects of this configuration were 
considered desirable. Although the inability of such a leading edge to 
support an attached shock is well known, the significance of its effect 
on impeller flow inducement is not yet fully recognized. Deviation from 
the conventional theoretical supersonic compressor inlet phenomena 
(ref. 3) appears more a difference of degree than character, since the 
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existence of subsonic axial-flow velocities before any type impeller 
precludes the exclusive existence on the blade trailing face of either 
compressions or expansions without subsequent waves of the opposite sign- 
in any case, however, the blade trailing— surface contour is probably 
still a contributing factor in the determination of the inflow velocity. 
One characteristic difference between the swept and unswept blade lies 
in the ability of the hub inlet flow to effect the rotor inlet flow at 
larger radii. No comprehensive experimental evaluation of the swept 
blade effects can be gathered from the current work because of the 
aforementioned flow restriction, which prevented the attainment of 
design inlet conditions. 

Flow limits. - A reasonably regular curve could be drawn through 
the maximum weight flow points of figure 6- The mininum measured weight 
flow points, however, do not allow a regular curve at all. The reason 
for this must be found in the inability to locate a surge point at any 
speed. Moreover, continual cognizance of the threat of structural 
f a j 1 nvp along with observation of the vibrations with increased back 
pressure prevented more elaborate exploration in the low weight flow 
regions. Only for the run at 1100 feet per second, where there was 
unusually smooth operation, is there any indication of a flow limit. 

The static -pressure distribution along the outer casing is shown in fig- 
ure 8 for maximum and minimum weight flow points at tip speeds of 500, 
1100, and 1200 feet per second. For the upper curve at 1100 feet per 
second which corresponds to a much lower flow than the minimum 1200 feet 
per second flow point, the ratio of the local static pressure to the 
inlet total pressure may be seen to be of the order of unity at the 
impeller tip leading edge indicating a generally stagnant region, 
although there was no indication of audible surge- The type of stall, 
•whether rotating or stationary, and its radial extent were not known. 
Because it was observed that any back pressure change which changed the 
outlet conditions is immediately felt upstream even near maximum flow, 
the impeller itself is clearly unchoked. 

The original test program for the transonic mixed-flow impeller was 
to include a complete stator investigation following the impeller tests- 
These plans were altered when design impeller performance could not be 
obtained. It was thought, however, that perhaps the installation of 
stators would make possible an improvement in impeller performance by 
turning the outlet flow in a direction which would increase the effec- 
tive flow area, thereby helping to remove the choked region in the exit 
annulus. Thus stators were designed and built for a comb ination of the 
theoretical and experimental rotor exit flow conditions. As indicated 
by the tailed symbols in figure 6, only scant reduction in the flow 
restriction was noted. This amount of improvement was felt to be dis- 
tinct but insufficient to justify a detailed stator program; hence no 
further reference will be made to the stators, which basically had the 
effect of a small reduction in back pressure on the impeller. 
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Comparison of design and experimental performance . - The inability 
to obtain design weight flow makes it impossible to determine conclu- 
sively whether the design procedure accurately predicts the experimental 
flow. No detailed analysis was made by which the exact angle of inci- 
dence could be determined. The design computations were made for con- 
ditions which would be compatible with essentially zero angle of inci- 
dence at design flow. 

For the current design the mean static -pressure distribution along 
a streamline \|r = constant varied from nearly constant at the tip to a 
rapidly falling static pressure along the hub boundary. Thus the energy 
addition at design speed would be added all in the form of kinetic energy 
if design flow was realized. This condition, which should lead to effi- 
cient rotor flow, may lead to severe stator starting problems because of 
the high Mach number and consequently high shock loss level (see, e.g., 
ref. 13) • Moreover, when these starting difficulties exist, even the 
advantage of a favorable mean gradient in the impeller is questionable. 
Therefore, inasmuch as the axial-flow transonic work of reference 6 has 
indicated that a moderate static -pressure increase still gives very high 
efficiencies in the transonic speed range, the use of an increasing mean 
static pressure everywhere through the impeller would provide better 
over -all performance. 

Reference to the incidence angles (table I) indicates that the 
measured values do not approach the design values (approximately zero) 
as the design speed is approached. Instead, the angles increase, indi- 
cating that the flow restriction becomes more severe with increasing 
speed, possibly because of the increasing boundary layer growth and 
losses through the impeller which allowed scant increase in weight flow 
for the increase in rotor output. 

Table I also indicates that M-^^ is equal to or greater than 
M]_ -f- at each minimum weight flow point except that for a tip speed of 

1200 feet per second, where data were obtained for only a small weight 
flow range. While it does not appear feasible to detail the procedure 
by which the rotor flow breaks down, the deterioration of the tip flow 
clearly precedes that of the hub, a phenomenon associated more closely 
with the stall or surge region than design operation, however. 

As shown in figure 9, the variation in both the theoretical and 
experimental total -pres sure ratios over the outlet annulus is large. 

This variation is probably detrimental from both the viewpoint of dif- 
ficulty of staging and efficiency of stator diffusion since the radial 
energy gradient would probably lead to excessive mixing in the diffusing 
process after the dynamic balance is relaxed. For this particular 
impeller, which was designed to investigate other aspects of impeller 
performance, the radial variation in total pressure was not considered 
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important. However, the required alteration to the design procedure 
could be made to maintain constant theoretical work input along the 
radius, particularly if the restriction to radial blade elements is eli- 
minated or relaxed. The total-pressure ratio or work output variation 
would then be a function of the efficiency variation along the blade 
span. 


Blade Structural Failure 

During a routine inspection after approximately 35 hours and 15 min- 
utes running time, one of the impeller blades was found to have a 1^ inch 
crack beginning at the leading edge and running just above and nearly 
parallel to the hub. The crack closely checked a nodal point located as 
shown by the inner sand pattern (fig- 4) during a preliminary check of 
the static vibrational characteristics of the impeller. The frequency 
of this mode of static vibration was observed to be between 1180 and 
1210 cycles per second. A small amount of undercut had previously been 
noted at the hub of each blade. This defect, coupled with the existence 
of several relatively easily excited vibrational modes, indicated the 
possibility of this type failure. 

The impeller was cut from a solid aluminum forging and hence pro- 
vided little internal damping. However, since the original crack plus 
a smaller one of the same type discovered later were small, it was felt 
that some further tests might be run if some form of damping could be 
installed. Consequently, a 0.030 inch hole was drilled in each blade in 
the tangential direction at a point midway between the hub and tip near 
the leading edge. Through the holes were strung 18 loops of 0.006 inch 
light filament wire. The two ends were cemented, one on each side of one 
of the blades (fig. 10) . This arrangement reduced the previously 
recorded static natural frequencies to approximately one -seventh their 
previous values. With this configuration the originally scheduled test 
points were completed in about 6 hours running time before the wire 
failed, probably because of the unequal tension on the loops. Little, 
if any, further deterioration of the blades was observed and no apparent 
aerodynamic effects of the wire could be detected. 


CONCLUDING REMARKS 

An analytical and experimental study of a high solidity transonic 
impeller was made. 

An approximate procedure for design of high solidity impellers 
assuming axially symmetric, nonviscous flow is developed and presented 
along with the design and experimental results for a 12 -inch-diameter , 
swept back leading edge, transonic, axial discharge impeller. 
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Test results at the design equivalent tip speed of 1200 feet per 
second indicated a maximum pressure ratio of only 2 . 164 at a weight flow 
of 26.47 pounds per second (33-71 lb/sec-sq ft) and an impeller adia- 
batic efficiency of 0.81. The measured performance was restricted par- 
ticularly at and near design speed to performance figures considerably 
below theoretical by a flow restriction in the outlet annulus. The low 
speed performance appeared very good. Partial structural failure of the 
impeller occurred before alleviation of the outlet annulus flow restric- 
tion could be accomplished. 

The effects of blade sweep back, though not fully understood, do 
not appear undesirable inasmuch as this impeller operated without audible 
surge in spite of very large incidence angles. 

Test results and analysis strongly indicate that for future designs 
a rise in mean static pressure through the impeller in order to reduce 
absolute discharge velocity and an approach to equal energy input along 
the rotor blade span are desirable. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX - SYMBOLS 

The following symbols are used in this report : 
velocity of sound, ft/sec 

specific heat at constant pressure, ft -lb/ (slug) (°P) 

impeller diameter, ft 
derivative with respect to time 
force per unit mass 
height of annular stream tube, ft 

absolute Mach number, ratio of absolute air velocity to local 
velocity of sound 

relative Mach number, ratio of air velocity relative to rotor to 
local velocity of sound 

streamline number 

number of blades 

total pressure, lb/sq ft 

stream pressure, lb/sq ft 

radius of curvature in radial-axial plane, ft 
compressor radius, ft 

streamline projection in radial-axial plane, ft 
stream temperature, °R 

velocity of rotor (car) at radius r, ft/sec 
absolute velocity, ft/sec 
velocity relative to rotor, ft/sec 
weight flow, lb/sec 

weight flow between adjacent streamlines 
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a incidence angle, angle between inlet flow and blade mean line at 

leading edge 

P angle between compressor axis and absolute flow direction, deg 

P' angle between compressor axis and air velocity relative to rotor, 

deg 

V ratio of specific heats 

S ratio of actual inlet total pressure to standard sea-level pres- 

sure P t ^q/2116 

^ad adiabatic efficiency 

9 angle between blade element at blade leading edge and subsequent 

element both projected on plane perpendicular to compressor 
axis, radians (see sketch preceding eq. (l)) 

'f& square root of ratio of actual inlet stagnation temperature to 

standard sea-level temperature , rj Tq/518.4 

p fluid density, slugs/cu ft 

T blade thickness 

cp upsweep angle, tan -1 V r /v z 

cp’ tan-1 V r /V zd 

t stream function along which w = constant 

cu angular velocity, radians/ sec 

Subscripts : 

df driving face of blade 

H hub 

m midpoint of a particular incremental annulus 

n points along a streamline 

r radial component 

rz radial-axial plane 
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T total 

t tip 

tf trailing face of blade 

z axial component 

zQ axial-tangential plane 

9 tangential component 

0 position at beginning of entrance section 

1 impeller inlet 

2 impeller exit 
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TABLE I. - PERFORMANCE DATA 


Equivalent 

tip 

speed, 

Ut 

Pressure 
ratio, U 

p 2 / p 0 

Adiabatic I 
ff iciency, 

^ad 

Equivalent 

weight 

flow, 

w<v/e 

6 

Tip 

angle 

of 

inci- 

dence, 

a t 

Hub) \l 
angle 
of 

inci- 

dence, 

a H 

Absolute P 
inlet 
Mach 
number 
at tip, 

M i,t 

ib solute 
inlet 
Mach 
number 
at hub, 

m 1,H 

1198 
1202 

1199 
1198 
1196 

2.16 
2.14 
2.10 
2.07 
2.06 1 

0.81 

.79 

.77 

.76 

.76 

26.47 

26.04 

25.19 

24.27 

23.42 

19°32' 
19°30 ' 
21°54 ' 
23°10' 
24°12 ' 

5°39 ' 1 
6°45' 
7°14 ' 
8°56 ' 
9°49 1 

0.54 

.54 

.48 

.45 

.43 

0.49 

.47 

.46 

.44 

.42 

1101 

1099 

1100 
1099 

2.00 
1.87 
1.79 
1.70 | 

0.86 

.76 

.71 

.64 

25.80 

22.83 

18.97 

14.73 

18°16 ' 
23°15 ' 
28°25' 
36°00' 

4°31' 
8°57' 
15°02 ' 
13°44 1 

0.53 
.41 
.31 
* . 17 

0.47 

.40 

.32 

.33 

958 

1000 

998 

999 
998 

1.79 

1.79 

1.68 

1.69 

1.63 

0.93 

.84 

.76 

.75 

.70 

24.85 

24.36 

22.94 

21.39 

19.57 

15°55 1 
17°53' 
20°57' 
23°41' 
26°18 1 

2°02' 
3°52' 
6°40' 
9°00' 
11°15 ' | 

0.503 

.480 

.42 

.37 

.32 

0.44. 

.43 

.39 

.36 

.33 

799 

799 

800 
799 
802 
799 
799 

1.51 

1.48 

1.49 
1.45 
1.47 
1.43 
1.41 

0.91 

.87 

.88 

.83 

.85 

.79 

.77 

23.59 

21.98 

20.93 

20.02 

18.91 

18.25 

16.37 

12°07' 
15°49 ' 
16°01 1 
19°33' 
18°44 ' 
23°53' 
25°32' 

-0°53' 
1°59' 
2°19' 
5°13' 
5°34' 
7°55' 
10°35 1 

0.49 

.42 

.41 

.36 

.37 

.29 

.27 

0.41 

.37 

.37 

.33 

.33 

.30 

.271 

700 

701 
700 
699 
699 
699 

1.39 

1.37 

1.35 

1.33 

1.31 

1.28 

0.95 

.90 

.87 

.81 

.77 

.70 

22.09 

20.20 

18.49 

16.65 

14.52 

12.40 

10°50 ' 
14°45 1 
17°44 1 
21°38 ' 
25°58 1 
28°47 ' 

-1°33' 
1°03' 
3°23 ' 
7°11' 
11°32 ' 
14°38 1 

0.44 

.38 

.34 

.28 

.23 

.19 

0.37 

.33 

.31 

.27 

.27 

.21 

500 

502 

500 

500 

499 

500 
500 
500 
500 
500 

1.19 

1.19 

1.21 

1.18 

1.19 

1.19 

1.20 
1.18 
1.18 
1.14 

1 0.95 

.96 
1.03 
.92 
.98 
.94 
1.00 
.92 
.88 
.72 

20.36 
19.72 

19.37 
17.46 
17.19 
15.53 
15.43 
14.25 
13.21 
10.31 

3°28' 
4°56 1 
3°45' 
9°14' 
7°50 
11°48 
13°34 
14° 51 
19°18 
25°13 

1 -7°37 ' 
' -6°43' 
' -7°45' 
' -2°26' 
' -4°22' 
' -1°21' 
' 0°35' 

' 1°40 

i 4°46 1 
' 9°53 

0.41 
.39 
.41 
.33 
' .35 

' .30 

' .28 

• .27 

1 .22 

» .17 

0.34 

.32 

.34 

.27 

.29 

.26 

.24 

.23 

.21 

.17 


f- Measured static pressure au jjuxuo 

An orifice l/2 inch upstream of the impeller inlet gave Mi,t 


0.17. 


Blade thickness, t/r , percent tip radius 
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Figure 1. - Blade thickness distribution along tip mean line. 


Blade relative velocity, V', ft/sec 
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Figure 3. - Velocity diagram for tip mean annul ax streamline. 
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Figure 4. - 12-inch transonic high' solidity impeller. 
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Figure 6. - Over-all performance of transonic mixed-flow 
impeller . 
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Figure 7. - Variation of hub and tip angle of incidence with impeller tip speed for 
minimum back pressure condition. 
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Figure 9. - Measured flow conditions at impeller 
exit for design speed and maximum weight flow. 
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Figure 10. - Transonic impeller after Installation of damping wires. 
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